Generalization of Child-Langmuir Law for Non-Zero Injection Velocities in a Planar 
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The Child-Langmuir law relates the voltage applied across a planar diode to the saturation value 
Jcl of current density that can be transmitted through it in case the injection velocity of electrons 
into the diode is zero. The Child-Langmuir current density Jcl is, at the same time, (i) the maxi- 
mum current density that can be transmitted through a planar diode, (ii) the current density below 
which the flow is steady and unidirectional in the long time limit and (iii) the average transmitted 
current density for any value of injected current density above Jcl- Existing generalizations of 
Child-Langmuir law to non-zero velocities of injection are based on the characteristics (i) and (ii) 
of Jcl- This paper generalizes the law to non-zero velocities of injection based on the characteristic 
(iii) by deriving an analytical expression for the saturation value of current density. The analytical 
expression for the saturation current density is found to be well supported by numerical computa- 
tions. A reason behind preferring the saturation property of the Child-Langmuir current density as 
the basis for its generalization is the importance of that property in numerical simulations of high 
current diode devices. 



I. INTRODUCTION 

The celebrated Child-Langmuir law is widely used in 
many areas of physics where charged particle transport 
across a potential difference is involved [1]. It gives the 
maximum current density Jcl that can be transported 
from an infinite planar cathode at zero potential to an 
infinite planar anode parallel to the cathode at a distance 
L and potential V when the initial velocity vq of electrons 
at the cathode is zero. This maximum value is a result of 
the space charge effect and, in the non-relativistic limit, 
is given by 



Jcl 



1 /2e\ 1/2 
L 2 9tt 



(1) 



Here, e is the magnitude of the charge on an electron and 
m is its rest mass. For convenience, the terms "current" 
and "current density" are used interchangeably. They 
refer to the current density denoted by J. 

There have been several attempts to generalize the law 
(1) to other situations such as to higher dimensions or to 
non-zero velocities of injection [2]- [6]. The characteris- 
tics of Jcl on which such generalizations are based are 
one or the other of the following: (i) Jcl is the maximum 
current that can be transmitted across one-dimensional 
planar diode with zero velocity of injection, and (ii) the 
flow of current in the long time limit is steady and uni- 
directional if the injected current is less than Jcl but 
oscillatory in time above it. The generalization based 
on characteristic (i) to non-zero injection velocity, Vo, of 
electrons in to a planar diode shows that the space charge 
limiting current 
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(2) 



is the maximum current that can be transmitted across 
the diode and is identified as the generalized Child- 



Langmuir current [5]. The generalization based on char- 
acteristic (ii), on the other hand, shows that the current 
in the diode is unidirectional till it attains the value Jbf 
and starts reflecting thereafter with 



Jbf = Jcl 
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This current, Jbf, marking the bifurcation of the state 
from completely to partially transmitting is proposed in 
Ref. [6] as a generalized Child-Langmuir current. 

A third characteristic that is commonly ignored in 
such attempts is the fact that Jcl is the saturation cur- 
rent that is transmitted for zero velocity of injection for 
any injected current higher than Jcl- It is the existence 
of saturation current in a planar diode that we wish to 
explore in this communication for non-zero velocities of 
injection and call the relation between such a saturation 
current, the injection velocity, the diode spacing and volt- 
age as the generalized Child-Langmuir law. 

The reason for emphasising the saturation property is 
that it helps in simulations involving high current de- 
vices. For, the existence of a saturation current implies 
that one need not bother about the material dependent 
cathode characteristics (the Richardson-Dushman law in 
case of thermionic emission or the Fowler-Nordheim law 
for field emission) so long as the current emitted from the 
cathode surface is more than what gives rise to the sat- 
uration value of the transmitted current. This condition 
simplifies considerably the task of numerical simulation. 

Though the injection velocity of electrons in many de- 
vices can be taken as zero, in several others non-zero in- 
jection velocities do arise (see [1] and references therein). 
The magnitude of the velocity can be related to physi- 
cally measurable quantities such as the temperature in 
case of thermionic emission or the frequency of the inci- 
dent field in case of photoemission. Since, as in the case 
of vq — 0, existence of a saturation current for vq =/= 



1 



may enable one to ignore the cathode characteristics for 
appropriate injected currents, it is important to know 
for what values of injected current and to what value, 
if any, does the transmitted current saturate if vo ^ 0. 
The main result of this paper is that as the current in- 
jected in to a planar diode is increased beyond its space 
charge limiting value, the current transmitted through it 
saturates asymptotically to the value 



Jgcl — Jcl 1 + 



mv. 
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(4) 



which we call the generalized Child-Langmuir current. 
The equation above is the desired generalization to v ^ 
of the Child-Langmuir relation (1) for v = 0. Note that 
Jgcl < Jbf < Jscl and that the equalities hold when 
v = 0. 

The relation (4) is derived in the framework of a phc- 
nomenological theory, called the classical theory. The 
validity of that theory is assessed by comparison with 
numerical solution of the equations. 

The paper is organised as follows. In Sec. 2 we recall 
the main results of the steady state theory leading to the 
expression for the space charge limiting current. Some 
details of derivation of those results are outlined in the 
Appendix A. The classical theory for currents exceeding 
the space charge limiting current is presented in Sec. 3. 
The results of numerical computation are presented in 
Sec.4. 



II. SPACE CHARGE LIMITING CURRENT 

Consider a diode consisting of two infinite parallel 
plates placed at z — and z — L where z is the direc- 
tion orthogonal to the plates. Let an electronic fluid of 
uniform and same cross-section as the plates move along 
the z-direction. Its motion in the Eulerian formalism is 
governed by the equations [7] 
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where p(z,t) and J(z,t) are respectively the charge and 
current densities of the electronic fluid, p(z, t) is the mo- 
mentum of the electron whereas <f>(z,t) and E(z,t) are 
the potential and the electric field at the position z at 
time t. 

We consider the case when a stream of electrons of 
constant charge density p enters the diode at z = 
with each electron in the stream entering with the same 
velocity vq along the z-direction so that 



p(0,t)=po, J{0,t) = p(0,t)v = J . 



(6) 



Furthermore, the plates at z = and z = L are assumed 
to be maintained at fixed potentials given by 



0(O,t)=O, <p(L,t) = V. 



(7) 



Because of the defining relation between the field and the 
potential, the condition (7) may be rewritten as 



E(z,t)dz = V. 



(8) 



A complete description of the motion of the electronic 
fluid is provided by solving equations (5) under the 
boundary condition (7) or (8). 

A state of interest, called a steady state, is the one 
that the system reaches asymptotically in time i.e. as 
t — > oo. The evaluation of space charge limiting current 
requires solving the steady state form of the equations 
obtained by equating to zero the partial derivatives with 
respect to time in the Eulerian formalism. It is evident 
from (5) that the current density in this case is not only 
independent of time but is also the same at every plane. 
If there is no reflection of electrons from anywhere in the 
diode then the current density in the diode is everywhere 
the same as is pumped in at its entrance so that, in the 
steady state, J(z, t) = p(z, t)v(z, t) = J . The first of the 
equations in (5) then assumes the form 



dz 2 



47rJo 



(9) 



The relation between <j> and v is provided by the steady 
state form of the law of conservation of energy, 
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where we have introduced the parameter e to denote the 
electronic energy and assumed that e is small enough for 
the motion to be considered non-rclativistic. On combin- 
ing (9) and (10) we arrive at the equation 



d 2 ^ 
dz 2 
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(11) 



where the scaled current a and the scaled potential <f> are 
defined by 



a = $ne\J \L 2 £ 

e 
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(12) 
(13) 



In terms of the scaled potential, the boundary condition 
(7) assumes the form 



(j)(z = 0,t) = 1, 4>{z = L,t) = V, 



where 



eV 

V = l + —. 

e 



(14) 



(15) 



Though the steady state may be determined by solving 
(11) along with the boundary conditions (14), it turns 
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out to be more convenient to derive it instead in the La- 
grangian formalism. We will, however, return to Eq.(ll) 
in the next section which concerns the situation involving 
returning currents. 

Alternative to the Eulerian formalism which charac- 
terizes the motion in terms of hydrodynamical variables 
like charge and current densities as function of space and 
time, is the Lagrangian formalism which is based on de- 
termining the trajectories of the electrons. An electronic 
trajectory in the one-dimensional motion under consid- 
eration is identified by the time t at which it enters the 
diode. In the non-relativistic limit (p — mv), the tra- 
jectory of an electron entering the diode at t = to is 
determined by the solution of the Llewellyn's equation 



d 3 z(t,t ) 

dt 3 



m 



dE (t) 
dt 



47rJ (t) 



(16) 



where E (t) = E(0, t) and Jo{t) = J(0, t) are the electric 
field and the current density at the entrance (z — 0) of 
the diode. The steady state solution of this equation 
has been derived in Ref. [5]. However, for the sake of 
completeness and clarity, and in order to fill the gaps 
in the available derivations, we present in Appendix A 
the solution of (16) for the steady state corresponding to 
unidirectional flow of current in terms of the following 
dimcnsionlcss variables 



z = 



E = 



V 
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t 
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The injected current Jo in the steady state is transmitted 
completely to the anode where the electrons arrive with 
velocity v\ which, by virtue of the energy conservation 
relation (10), is given by 



= Vv. 



(18) 



The main results of the steady state solution are: 



1 . The steady state is determined by solving the cubic 
(A9) for the transit time T(l) of an electron across 
the diode. The flow of current is unidirectional and 
independent of time if that cubic admits physically 
acceptable solution. 

2. The cubic admits two positive roots, 7^(1) and 
T 3 (l), given in (A14), if the value of the injected 
current a < ascL where 



ascL 



(i + Vv) 3 



(19) 



is the space charge limiting current. By invoking 
(12), it may be verified that the expression for the 
injected current Jo corresponding to a given above 
is the same as in (2) and that if v = then 



Jscl = Jcl, 



v Q = 0, 



(20) 



where Jcl is given by (1). In other words, for zero 
velocity of injection, the space charge limiting cur- 
rent is the same as the Child-Langmuir current. 



3. For < a < odr, the root 73(1) corresponds to 
negative velocities in certain region of the diode 
and hence is unacceptable in the said range of the 
injected current. However, both the roots, 7^(1) 
and T 3 (l) are acceptable if «dr < a < «scl where 



a UR = (i + v 3 /y 



(21) 



The suffix DR on a stands for " double root" : it is 
the value of the current above which steady state 
equation has two admissible roots. The linear sta- 
bility analysis shows that, in the range mentioned 
above, 7^(1) leads to unstable steady state. Hence, 
7^(1) defines the acceptable stable steady state for 
all a < a SCL . 

4. The flow of current ceases to be unidirectional and 
time-independent if the injected current exceeds 

The above-mentioned characteristics of the steady state 
solution are exhibited in Fig.l by the plot of the elec- 
tric field E at the cathode evaluated using (A16) for 
V = 1.2. The point A in the figure corresponds to the 
space charge limiting current ascL whereas the point B 
represents «dr- The branches OA and AB in that fig- 
ure correspond, respectively, to 72(1) and 7jj(l) given by 
(A14). The branch AB represents unstable steady state. 
The uppermost branch of the curve for a > «dr is a 
result of the classical theory outlined in the next section. 
The value of E at ckdr where the unstable steady state 
corresponding to 73(1) merges with the classical theory, 
found using (A16), is given by 



E BR= 3(1 + ^ 3 / 4 ). 



(22) 



We will see that, as it should be, this value is the same 
as that obtained by the classical theory. 

By transforming back to the unnormalizcd variables, 
the results above may be applied to the case of u = 0. 
It may be verified that the expression for space charge 
limiting current density Jo, found using (19) in this case 
is the well-known Child-Langmuir current Jcl given in 
(1) and that, correspondingly, E — 0. We thus see that, 
contrary to the standard treatments of the case of vo = 0, 
Eo = is not an imposed condition but is a natural 
consequence of the governing equations of motion. This 
fact has also been emphasised in [5]. 

Since the space charge limiting current oscl is the 
maximum current that can be transmitted across the 
diode, it is proposed in [5] that q;scl may be identi- 
fied as the generalized Child-Langmuir current and (19) 
as the generalized Child-Langmuir current-voltage rela- 
tion. However, as we have argued in the Introduction, 
an important characteristic of Child-Langmuir current is 
that it is the current transmitted across the diode for 
any injected current above «scl- We, therefore, need 
to examine the behaviour of the transmitted current for 
a > a S cL- 

As shown in the Appendix A, the system does not ad- 
mit a steady state corresponding to unidirectional flow of 
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current if a > «scl- Hence, the behaviour of the system 
for a > ascL is determined by solving time-dependent 
equations. Since solving time-dependent equations ana- 
lytically is a formidable task, those equations are solved 
numerically. However, a qualitative understanding of the 
phenomenon is achieved by a simple model which as- 
sumes that, though for a > q;scl the current may not 
be unidirectional, it nevertheless attains a steady state 
with current flowing in two directions in some region in 
the diode. This theory, called the classical theory, is out- 
lined in the next section. It is along the lines of the theory 
in Ref. [7] for the case of grounded plates. 



III. ASYMPTOTICALLY SATURATING 
CURRENT 

The analysis outlined in the Appendix A holds good 
for a < ascL m which case the injected and the trans- 
mitted currents are equal. The classical theory assumes 
that, beyond ascL, the minima of the potential within 
the diode assumes a value which is such that the velocity 
of an electron there reduces to zero. The equations (10) 
and (13) show that v = is attained when (f> — 0. Hence, 
the boundary condition at the position z — z m of the 
potential minima is 



dz 



(z = z m ). 



(23) 



After the electrons come to rest at z — z m , it is assumed 
that a fraction / of them is transmitted while the remain- 
ing fraction 1 — / is reflected back towards the cathode. 
Thus, in the region < z < z m , there are two currents: 
the injected current a moving away from the cathode and 
(1 — f)a moving towards it. In the region z m < z < 1, 
there is unidirectional flow of current of magnitude fa 
towards the anode. Note that the net current at every 
plane in the two regions is fa. Since the charge densi- 
ties at a point due to different currents at it add, the 
Poisson's equation (11) in the two regions assumes the 
forms 



d 2 ^ L = 4a(2 - f) 
dz 2 



d 2 c 




dz 2 



< z < z r . 



Zm < Z < 1. 



Both these equations are of the form d 2 <j)/dz 2 
and can be recast as 



_d 

dz 



dz 



dz 



(24) 
(25) 
C/^ 2 

(26) 



The solution of this equation, satisfying the boundary 
condition (23) reads, 



dz 



(27) 



Now, recall that in a one-dimensional electrostatic situa- 
tion involving charges of only one sign, there can be only 



one extrema. In the problem at hand, the minima at 
z — z m , therefore, is the only one admissible so that the 
electrons are decelerated while moving from the cathode 
to the position of minima at z = z m and are accelerated 
thereafter. Hence, the minus sign in the solution above 
holds in the region < z < z m whereas the plus sign is 
to be retained in the region z m < z < 1. On integrat- 
ing (27) in the two regions with appropriate values of C 
along with the boundary conditions at z = 0, 1 given in 
(14), it follows that 



fl* = 1 - zy/a{2-f), 
fl' = {z-\)Jal+V*l\ : 



< Z < Z m , 
m<Z<l. 



(28) 
(29) 



The electric field at the cathode, obtained by substituting 
(28) corresponding to z = in (27) reads 



V"(2 - 7) 



(30) 



The unknowns z m and / in the expressions above are de- 
termined by demanding that cf> in the two regions should 
match at z = z m . To that end, note that, by virtue of 
(23), the expressions (28) and (29), evaluated at z = z m , 
lead to the relations 



1 



^3/4 

7^7' 



(31) 
(32) 



These equations determine two unknowns, the position 
z m of the plane of reflection and the fraction / of the 
injected current that is transmitted. On eliminating z m 
between these equations, the expression determining / 
reads 



\[a = 



1 



1/3/4 



(33) 



It is straightforward to verify that d//da < 0. Hence, / 
is a decreasing function of a. Now, (33) shows that the 
value / = 1 is attained when a = «bf where 



«BF= (l + ^ 3/4 ) 2 . 



(34) 



According to classical theory, this is, therefore, the mini- 
mum value of a for which reflecting solution exists. The 
subscript BF on cubf indicates that the steady state bi- 
furcates into two different types at q;bf : the steady state 
corresponds to unidirectional flow of injected current a if 
a < aBF whereas it starts reflecting if a > ckbf • On com- 
paring (34) with (21) it follows that «dr = «scl i.e. the 
bifurcation of the steady state from unidirectional flow to 
partially reflecting flow takes place at the same value of 
current above which there are two solutions for the uni- 
directional flow, one stable and another unstable. The 
electric field at the entrance of the diode for a = cubf, 
obtained using (30), reads 



4 



E, 



OBF 



dz 



= (1 + ^3/4) (35) 
-0 o 



This is the same as the one in (22) obtained using the 
steady state solution for a = a DR , 

It may be verified that the current Jo corresponding 
to Q!bf is given by (3). Since the bifurcation point marks 
the onset of reflections, which is one of the character- 
istic properties of the Child-Langmuir current, the cur- 
rent Q!bf at the bifurcation point is identified in Ref. [6] 
as the Child-Langmuir current. However, the numeri- 
cal solution of equations shows that reflection does not 
start at the bifurcation point. It starts only when the 
injected current exceeds its space charge limiting value 
«SCL which is higher than its value at the bifurcation 
point except when v = in which case, a comparison of 
(2) and (3) shows that 



Jbf = Jscl, v — 0. 



(36) 



On comparing with (1) it follows that these currents are 
the same as the Child-Langmuir current Jcl- 

The behaviour of Eq for partially reflecting solution 
is exhibited in Fig.l by the branch of the curve above 
the bifurcation point B. That branch is obtained by 
computing first the fraction / of the transmitted current 
by solving 33 for a given a. The value of / so obtained 
is substituted in 30 to evaluate E . 

Before proceeding further, we digress to rewrite 
Eq.(32) to express the transmitted current fa in the form 



fa = 



(1 + eV /e) 
(1- 



3/2 



(37) 



This resembles the expression derived by Langmuir [8] for 
the thermionic current from a hot cathode emitting into a 
diode. Langmuir's expression in question reduces exactly 
to (37) if the cathode temperature is zero whereas it gives 
the current approximately if the temperature is non-zero 
but low. It should, however, be emphasised that in the 
model under discussion, all the electrons are injected in 
to the diode with the same velocity whereas, in the pa- 
per of Langmuir, the electronic velocities at the entrance 
are distributed according to the Boltzmann distribution. 
The position z m of the potential minimum in (37) deter- 
mining fa itself depends on /. Hence, evaluation of (37) 
requires a relation between z m and / which, in our case, 
is contained in (31). 

The Eq.(33) leads to the following useful expression 
for the transmitted current, 



«tr = fa = 



y3/4 + 



2-/ 



(38) 



for all a > «bf- In terms of the unnormalized variables, 
this reads 



Jtr = Jcl 



2\ 3 / 4 
WIVq \ 

2eVJ 



f 



+ 1 



mv. 



3/4' 



2eV J 



for all Jo > Jbf- Since < / < 1, this expression, in the 
limit of vanishing injection velocity (e — > 0) yields 



Jtr = Jcl, 



vq = 0, Jo > Jbf, 



(40) 



where Jcl is the Child-Langmuir current (1). This is 
the current transmitted for any injected current above 
its space charge limiting value. This value of current is 
the same as the value of the space charge limiting current 
(2) corresponding to vq = 0. It therefore follows that, for 
vo = 0, the transmitted current saturates at Jcl as soon 
as the injected current attains its space charge limiting 
value. 

The saturation current for non-zero injection veloci- 
ties, obtained from (38) in the limit f —* 0, reads 



«tr — ► accL, 



as / -> 0, 



(41) 



where accL is the saturation current given by 

a G CL - V z ' 2 . (42) 

We call the saturation current in the equation above as 
the generalized Child-Langmuir current. Its unnormal- 
ized form is given by (4). 

The finite value of the transmitted current fa in the 
limit / — > implies that it is achieved when a — > oo. 
Thus, the current transmitted across the diode saturates 
at the value given by (42) when the current injected in it 
is increased indefinitely. 

We assess the validity of the classical theory by com- 
paring its predictions with numerical results in the next 
section. 



IV. NUMERICAL RESULTS 

In this section we present a comparison of the predic- 
tions of the steady state theories and numerical compu- 
tations carried out using the sheet model [7]. The com- 
parison is carried by computing the transmitted current 
Jtr as a function of the injected current Jo (a is related 
to J by Eq. (12) ). 

The injected current Jo in the sheet model corresponds 
to injecting TV sheets per unit time (L/(2eV/m) 1 ^ 2 ) 
into the diode, each having a surface charge density 
a = Jq/N. By counting the sheets starting from the 
one closest to the collecting plate, it can be shown that 
when there are M sheets in the drift space, the electric 
field on the k th sheet at the position Zk from the entrance 
plate is given by 



E(z k ) = -Atto 



1 M 1 



i=i 



+ 



V 



(43) 



The (unnormalized) potential at any point z between the 
(k) th and (k + l) th sheet is similarly given by 



(39) </>{z) 



Vz AttJq 
1 N~ 



/ M \ M 

kL-J2 z i)+ L Yl Zi 

\ 1=1 / l=k+l 



(44) 
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where the sheet number is counted from the anode (right) 
plate. The force equation for the sheet is thus 

m ^ = eE{Zk) (45) 

The numerical procedure consists of solving the coupled 
set of equations for Zi(t), (i = 1,2,.. .) as a function of 
time by choosing the time-step At of integration and the 
rate of injection N of the sheets in such a way that desired 
accuracy is achieved. Typically, At — 8.532 x 10~ 13 s, 
N = 1000 for J = 20000 A/m 2 ' 

The numerical computations confirm the theoretical 
prediction that, for a < ascL, the system approaches a 
steady state asymptotically as t — * oo. However, contrary 
to the premise of the classical theory, the asymptotic be- 
haviour of the system in time is found to be oscillatory 
if a > c^scl- Hence, for a > ascL, the predictions of the 
classical theory are compared with the results of numer- 
ical computations averaged over several cycles of oscilla- 
tions after the initial transients have settled. See Ref. [9] 
for a plot of oscillating transmitted current as a function 
of time for injected current in excess of the SCL current 
in a short-circuited diode. 

The Fig. 2 is a plot of theoretical predictions and nu- 
merically computed values of the transmitted current 
Jtr as a function of the injected current Jo for e = 
and for e/eV = 0.2 with L = .008m and V = 250V. 
The bold straight line represents Jtr = Jo which is the 
prediction of the steady state theory for currents below 
the space charge limiting value. The square dots on that 
line are the results of numerical computations indicating 
agreement between analytical and numerical solutions for 
currents below the space charge limiting value. 

The numerical results of computation for currents ex- 
ceeding the space charge, carried in the manner outlined 
above, are represented by the square dots. The dashed- 
dotted straight line (a) is drawn to join the dots corre- 
sponding to e = whereas the dashed-dotted curve (b) 
for e/eV = 0.2 joins the dots corresponding to e = 0.2. 
Recall that the classical theory predicts that if e = 
then Jtr = Jcl for all Jo. For the values of the param- 
eters used in the figure, the value of the Child-Langmuir 
current is Jcl — 143.9 which practically coincides with 
the corresponding numerical results represented by the 
straight line (a). Hence, the theoretical predictions are 
in agreement with numerical results for any value of the 
injected current if e = 0. 

The theoretical plot for e/eV = 0.2 is represented by 
the dashed curve marked "classical theory". It is ob- 
tained by solving (33) for / for a given value of a and 
using the value of / so obtained to evaluate the trans- 
mitted current Jtr = fa- On comparing the classical 
theory plot with the corresponding results of numerical 
computations represented by the curve (b), it is seen that 
the agreement between the two is not close. For, in dis- 
agreement with the monotonic behaviour predicted by 
the classical theory, the numerical plot exhibits a max- 
ima in cutr as a function of a. The reason behind the 
appearance of the maxima is left as an open question. 
However, the asymptotic limit of the numerical plot of 



cktr appears to be in agreement with the prediction (41) 
of the classical theory. In order to confirm the apparent 
agreement we examine the asymptotic behaviour of the 
curve 

Jtr = a + bJ~' 3 (46) 

fitting best the numerical results. The Fig. 3 shows some 
values of Jtr computed numerically along with the curve 
(46) fitting best the numerical results as a function of 
J . The best fit is obtained for a = 193.58, b = 2318.83 
and (3 = 0.40. Hence, (46) approaches asymptotically the 
value Jtr = 193.58 which agrees very well with the value 
Jtr = 189.37 obtained using the theoretical expression 
(42). We have observed similar agreement between the- 
oretical prediction and the numerical computations for 
several other values of the injected velocity. 

Finally, we study the behaviour of the (time-averaged) 
position of the virtual cathode as a function of the in- 
jected current J for e = 0.2. Fig. 4 shows a plot of 
the time-averaged normalized position of the potential 
minimum (z m ) as a function of J . As in the case of 
short-circuited diode, z m decreases with J . Fig. 4 also 
shows the prediction of the classical theory (bold line) 
and the best fit of the function ai + biJ^ 5 (dashed line) 
for J in the interval [1000,20000]. The value of the pa- 
rameters are : a\ = —0.007, b\ = 1.384 and 5 = 0.335. 
Clearly, the position of the virtual cathode moves to the 
injection plane as Jo — > oo. 

V. CONCLUSIONS 

In this paper, we have sought to generalize the Child- 
Langmuir law for non-zero injection velocities based on 
the saturation property of the time-averaged transmitted 
current. The reason for preferring the saturation prop- 
erty for the said generalization is its importance in nu- 
merical simulations. For, as stated in the Introduction, 
the emission characteristics of cathode and the details of 
the emission mechanism do not play any role in determin- 
ing the transmitted current if that current is independent 
of (or weakly dependent on) the injected current. This 
is the case when the current emitted from the cathode 
surface (J ) is more than what gives rise to the satu- 
ration value of the transmitted current. This condition 
simplifies considerably the task of numerical simulation. 

Thus, if we characterize the Child-Langmuir current 
as the saturation current then ckgcl given in Eq. (42) is 
the Child-Langmuir current generalized to non-zero ve- 
locities of injection. That expression for the saturation 
current is derived by assuming a steady flow of current 
under the condition that a part of the current is reflected 
from the virtual cathode. Eq. (42) is found to be in 
fairly good agreement with numerical computations. In 
contrast to the case of zero velocity of injection for which 
saturation is attained at the space charge limiting value 
of the injected current, the saturation in the case of non- 
zero velocities is attained only asymptotically. Hence, 
though the details of the emission characteristic of cath- 
ode may be ignored safely in numerical simulations for 
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any injected current above its space charge limiting value 
in case vq = 0, the injected current for vq =^ needs to 
be very much above its space charge limiting value for 
cathode characteristics to become unimportant. 



APPENDIX A: STEADY STATE IN 
LAGRANGIAN FORMALISM 



Now, if T(l) is the time that an electron takes to reach 
the collector plate after its entry in to the diode then the 
Equations (A5) and (A6) yield 



vi = l-^f(l)£b + ^ 2 (1), 
1 = T(1)-V(1)£„ + ^ 3 (1)- 



(A7) 
(A8) 



In this Appendix we present steady state solution of 
the Eq.(16) describing the flow of current in Lagrangian 
formalism. 

Integration of (16) yields 



dv(t,t ) 
dt 



e 
m 



E (t)+4n [ drJo(r) 

Jt n 



(Al) 



This equation is arrived at by noting that the electron 
at time t = to is at z — 0. Hence, the initial condition 
is the same as the condition at z = so that, by virtue 
of the last equation in (5), dv(t,to) / dt = —eEo/m at 
t = to. The equation (Al) may be integrated further 
to obtain expressions for v{t, to) and z(t,to). Like (Al), 
those equations involve two unknown functions: the field 
and the current at all times at the entrance plate. The 
current at the entrance plate consists of the current being 
constantly pumped in and the current returning, if any, 
after reflection from some position between the plates. 
Our interest here is in the steady state without returning 
currents. The current at the entrance plate in that case 
is the same as that being pumped in i.e. 



Jo(t) = Jo = povo = -\po\vo- 



(A2) 



We normalize the variables as in (17) and integrate (Al) 
in the absence of returning currents to get 



If a - 

v(i, to) = 1 - - j_ dr E (T)dr + y(i-io) 2 , 



z(t,t ) = (t-to) 



1 

2./, 



dr {t-T)Eo(r) 



(A3) 



(A4) 



The unknown function Eo(t) in these equations is found 
by imposing the boundary conditions (7). That task is 
simplified in the steady state because the field E (t) at 
the entrance plate is then independent of time reducing 
(A3) and (A4) to the forms which are functions only of 
the time difference t — t : 

v(t,t ) = 1 - l(t - t )E + ^(t-io)\ (A5) 

z(t,i ) = (t - to) - tofEo + ^{t-iof. (A6) 

Note that these equations already incorporate the bound- 
ary condition at z = 0. The boundary condition (7) at 
z = L, restated in terms of the velocity by invoking the 
conservation law (10), leads to the expression (18) for 
velocity at the collector plate in terms of its potential. 



The two equations above determine the unknowns Eq 
and T(l). On eliminating Eq between these equation, 
the transit time T (1) is found to obey the cubic 

T 3 (l)- ^(i + o 1 )f(l) + ^ = 0. (A9) 

Of course, only the real positive roots of this equation are 
physically acceptable. The electric field at the entrance 
plate, obtained by substituting an acceptable solution of 
(A9) in (A7) is 



Eq 



T(l) 



a 



r T 2 (l) + l-^ 



(A10) 



Now, the nature of the roots is determined by the sign of 
the discriminant 



where 



so that 



27 

P = (1 + vi) 

a 



1 

— I 

27 



q = 



54 



a 



A = 



(AH) 



(A12) 



(A13) 



If A < then the three roots are real whereas one root is 
real and the other two complex conjugate of each other 
if A > 0. 

Note that A < and hence the three roots are real if 
a < (1 + vi) 3 . The roots in this case are 



ri(i) 


= -2A 1 / 3 cos 


r 2 (i) 


= 2A 1 /3 CO s( 


r 3 (i) 


= 2A 1/3 cos ( 



(A14) 



where 



A = 27 



-(l + «i) 

a 



3/2 



lau(0) = ( i(l + wi) 3 -l 



o<e<-. 



(A15) 



The root Ti(l), being negative, is physically unaccept- 
able. Hence, T 2 (l) and T 3 (l) are the admissible times of 



7 



transit. The expression for the electric field (A10) at the 
entrance of the cathode due to the positive roots assumes 
the form 



En 



3(1 + v^T) 



cos 



4(1 + vi) cos 2 



7T±£ 

3 

7T± 



'"1 



(A16) 



We now examine the question of admissibility of T 2 (l) 
and 23(1) based on other physics considerations. We 
will see that the slower of the two transit times, 7^(1), 
corresponds to negative velocity in some region if < a < 
(1 + v^ 2 ) 2 and to positive velocity all through the diode 

but an unstable steady state if (1+v^ 2 ) 2 < a < (1+vi) 3 . 

To that end, let v m be the minimum velocity at some 
place in the diode and let T m be the time that an electron 
takes to reach that position so that dv/dT = at T — 
f m . On solving (A3) (with i- t = f) for dv/df = 0, it 
follows that 



9 - 
T m = -7—E a , 

4a 

v m = i- JLisg. 



(A17) 
(A18) 



This equation shows that the minimum velocity is neg- 
ative if either E < — 4 v / a/3 or E > 4y/a/3. However, 
(A5) shows that if E < then v > 0. Hence, the condi- 
tion for the velocity to become negative in some region 
in the diode is E > 4 v / a/3 which, in conjunction with 
(A10), reads 



^T(i) - 1 - vsr) (^(i) - 1 + v^r) > 0. (aw) 



The positivity of T(l) implies that (A19) can be satisfied 
if 



f(l) > — (1 + v^T). 
\/ot 



(A20) 



This condition, for T(l) = T 2 (1),T 3 (1) given in (A14)), 
assumes the form 



cos 



7T± 



-) >\(l + V^){l + v 1 )- 1 ' 2 . (A21) 



It is easy to check that this implies cos((7r ± 9)/3)) > 
1/2 which, in turn, implies (ir ± 9)/3 < ir/3. Since 9 
is positive, the last inequality can be satisfied only by 
7r — 9. Thus, the minimum velocity can become negative 
only for 7^(1). Hence, 7^(1) is the only acceptable root 
for those values of a for which T 3 (l) leads to negative 
velocities. Those values of a may be identified by noting 
that the condition (A21) for T 3 (l) implies 



7T - 9 < 3 cos 1 



(1 + + 



-1/2 



(A22) 



Using the trigonometric relation cos(3x) = 4 cos 3 (a;) — 
3cos(a;), and the defining relation (A15) of 9, it is 
straightforward to show that (A22) holds if 



. ri . -3/2,2 

a < (1 + v{ ) = a D R, 



(A23) 



where «dr is given by (21). In other words, if a satisfies 
the condition (A23) then the root 7^(1) satisfies the con- 
dition (A20) required for the velocity to become negative 
in some region in the diode leaving 7^(1) as the only ac- 
ceptable root in the range < a < «dr whereas both, 
T 2 (l) and 73(1), are acceptable if «dr < a < (1 + vi) 3 . 
The ttDR denotes the value of the current above which 
there are two roots. However, the linear stability anal- 
ysis shows that the root 73(1) is unstable. Hence, only 
T 2 (l) corresponds to a stable steady state in the range 
< a < (1 + vi) 3 . 

The discriminant A>0ifa> (l+^i) 3 . The cubic 
(A9) determining the transit time in this case has one 
real and a pair of complex conjugate roots. Since the 
transit time is a real quantity, the real root is the only 
one which needs to be examined further for its accept- 
ability. To that end, we note that the product of the 
roots of (A9), given by — 54/a, is negative. Hence, the 
real root in case of A > must be negative. Since the 
time of transit is a positive number, the negative root is 
physically unacceptable. The cubic (A9), therefore, does 
not admit physically acceptable solution i.e. the system 
does not approach a time-independent steady state hav- 
ing unidirectional flow of current if a > q;scl where 



"SCL = (1 + Vi) S 



(A24) 



is the space charge limiting current: The space charge 
limiting current is the value of the current above which 
the flow of current is no longer unidirectional and time- 
independent. It is also the maximum current that can 
be transmitted across the diode even when the injected 
current is increased indefinitely. 

The solution of (A9) for a = q;scl gives the transit 
time as 



7scl(1) = 



which, on substitution in (A10) yields 



(A25) 



(A26) 



This is the electric field at the entrance of the diode when 
the value of the injected current a — ascL- 
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FIG. 1. The electric field at the cathode Eo plotted as a 
function of the current a for e/eV = 0.2. The branches OA 
(stable) and AB (unstable) correspond respectively to T2(l) 
and 73(1) given by (A14). The uppermost branch of the curve 
for a > ckdr(— Qbf) (point B) is a result of the classical 
theory. 
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FIG. 3. A fit of the transmitted current density using 
Jtr — a + bjq 13 as a function of the injected current den- 
sity Jo. The current density is units of Amp/m 2 . The solid 
line is the best fit while the squares mark the numerically 
computed values of the transmitted current using the sheet 
model. 
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FIG. 4. A plot of the time-averaged normalized position of 
the virtual cathode (solid squares) as a function of Jo- Also 
shown are the prediction of the classical theory (bold line) 
and the best fit using the function 01 + 61 J^ s (dashed line). 
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FIG. 2. The transmitted current density, Jtr, plotted 
against the injected current density, J , for two different value 
of e/eV with L = 0.008m V = 250V. The current density is in 
units of Amp/m 2 . The square dots are the results of the nu- 
merical computation. The bold line is Jtr = Jo- The straight 
line (a) is drawn to join the results of numerical computations 
for e = while the curve (b) joins the numerical results for 
e/eV = 0.2. The dotted line marked "classical theory" is the 
prediction of classical theory for e/eV = 0.2. 
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